X-ray observations of a new unusual magnetar Swift J1834. 9— 0846 



Oleg Kargaltsev 1 , Chryssa Kouveliotou 2 , George G. Pavlov 3,4 , Ersin G6gii§ 5 , Lin Lin 5 , 
Stefanie Wachter 6 , Roger L. Griffith 6 , Yuki Kaneko 5 , George Younes 2,7 

ABSTRACT 

We present X-ray observations of the new transient magnetar Swift J1834.9- 
0846, discovered with Swift BAT on 2011 August 7. The data were obtained 
with Swift, RXTE, CXO, and XMM-Newton both before and after the outburst. 
Timing analysis reveals single peak pulsations with a period of 2.4823 s and an 
unusually high pulsed fraction, 85% ± 10%. Using the RXTE and CXO data, 
we estimated the period derivative, P = 8 x 10~ 12 s s -1 , and confirmed the high 
magnetic field of the source, B — 1.4 x 10 14 G. The decay of the persistent X- 
ray flux, spanning 48 days, is consistent with a power law, F oc t -0 ' 5 . In the 
CXO/ACIS image, we find that the highly absorbed point source is surrounded 
by extended emission, which most likely is a dust scattering halo. Swift J1834.9- 
0846 is located near the center of the radio supernova remnant W41 and TeV 
source HESS J1834-087. An association with W41 would imply a source distance 
of about 4 kpc; however, any relation to the HESS source remains unclear, given 
the presence of several other candidate counterparts for the latter source in the 
field. Our search for an IR counterpart of Swift J1834. 9-0846 revealed no source 
down to K s ~ 19.5 within the 0V6 CXO error circle. 

Subject headings: X-rays: individual (Swift J1834. 9-0846) — gamma-rays: indi- 
vidual (HESS J1834-087) — ISM: individual (W41) — stars: neutron — X-rays: 
ISM 
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Introduction 



The population of magnetars has been growing rapidly in the last five years, reach- 
ing 24 objects as of August 2011. Originally comprised of Soft Gamma Repeaters (SGRs) 



and Anomalous X-ray Pulsars (AXPs) (Woods & Thompson 2006), the magnetar popula- 



tion now includes a few more neutron star (NS) groups that have been acknowledged as 
magnetar candidates. Most of these NSs are slow rotators emitting multiple, very short (a 
few times 100 ms) hard X-ray/soft 7-ray bursts. Their X-ray luminosities are likely pow- 
ered by the decay of their high magnetic fields (up to B ~ 10 15 G), rather than rotational 



energy losses due to their gradual spin-down (Paczyhski 1992 Duncan & Thompson 1992 



Thompson & Duncan||1995 1996). The current synergy between NASA's three observatories 



(RXTE, Swift, and Fermi) has enabled a much higher rate of discovery of these objects in 
the last three years. During July - August 2011 alone, two new candidate magnetars were 
discovered in X-rays, Swift J1822. 3— 1606 and Swift J1834. 9— 0846, when they triggered the 
Swift/Burst Alert Telescope (BAT) and the Fermi/Gamma-ray Burst Monitor (GBM). Their 
timing properties were subsequently established with Rossi X-ray Timing Explorer (RXTE) 
observations, clinching their magnetar nature. We report here on the X-ray spectral and 
temporal properties of the latter source. 

Swift J1834.9— 0846 was discovered on 2011 August 7, when a soft, short burst from the 
source triggered the BAT at 19:57:46 UT ( D'Elia et al.||20lT ); approximately 3.3 hours later, 
at 23:16:24.91 UT, another SGR-like burst triggered GBM from the general direction of the 
earlier BAT location ( Guiriec et al.||2011 ). Although the GBM location included a large area 
with several magnetar sources, the near time coincidence and the X-ray properties of these 



events pointed to a common origin of a new source (Barthelmy et al. 2011). The source 



triggered the BAT again on 2011 August 30 at 23:41:12 UT (Hoversten et al. 2011). 



Optical observations of the field ~ 16 min after the BAT trigger with the Special Astro- 
physical Observatory (SAO)/Big Telescope Alt-azimuth (BTA) 6-m telescope detected an 
object at magnitude R c = 23.44 ± 0.34 (Moskvitin et al. 2011). Simultaneous observations 



with the 1.5-m Observatorio de Sierra Nevada (OSN) telescope in the I band did not detect 
that object to a limit of / = 21.6 (Tello et al. 2011). Archival IR images of the region as 



part of the UKIDSS Galactic Plane Survey (Lucas et al. 2008) in the J, H and K bands 



on 2007 May 10, revealed two sources close to the Swift/X-ray Telescope (XRT) location of 



Swift J1834. 9— 0846 (Levan & Tanvir 2011). None of these objects coincided with the very 



precise X-ray position subsequently derived from our Chandra Target of Opportunity (ToO) 



observation (G6gu§ et al. 2011b). 



RXTE/PCA observations of the source on 2011 August 9-10 detected a coherent pul- 
sation at v = 0.402853(2) Hz, which corresponded to a spin period P = 2.482295 s (G6gii§ 
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on 2011 August 22 ( 


36gu§ et al. 


2011b) 



a time span of two weeks revealed a spin-down rate i> = —1.3(2) x 10 Hz s (Kuiper & 



Hermsen 2011). The corresponding estimate of the surface magnetic field, B = 1.4 x 10 14 



G, confirmed the magnetar nature of Swift J1834. 9— 0846. 

Swift J1834.9— 0846 is located in a field rich with high-energy sources, which include 
SNR W41 ( |Shaver fc Gosslp70| |Tian et al.|[2007| , the TeV source HESS J1834-087 QAharo- 



nian et al.|2005| ), the GeV source 2FGL J1834.3-0848 ( |Abdo et al.|20lT| , and the PSR/PWN 



candidate XMMU J183435.3-84443/CXOU J183434.9-084443 QMukherjee et al.|[2009| [MiF 



anovic et al. 2011). Attempts to understand the nature and relations between these sources 
had already prompted X-ray observations with CXO and XMM-Newton before the discovery 
of Swift J1834.9-0846 flMukherjee et al.||2009[ [Misanovic et al.|2011[ ). We have triggered ad- 
ditional observations of the region with both CXO and XMM-Newton. Here we describe the 
analyses of the RXTE, Swift, Fermi, and CXO data and compare them to the earlier obser- 
vations. The XMM-Newton results will be reported in a separate paper. Section [2] describes 
the data sets presented here, and Section [3] presents the CXO location and discusses possible 
optical counterparts. We present the lightcurve of the persistent emission in Section [4] and 
the results of our timing and spectral analyses in Sections [5] and |6j respectively. Finally, we 
compare the properties of Swift J1834.9— 0846 with those of other magnetars and discuss the 
possible relation of Swift J1834. 9— 0846 to other sources in the field in Section [7j 



2. X-ray observations and data reduction. 

The field of Swift J1834. 9— 0846 was observed in X-rays on 29 occasions with several 
telescopes; the majority was in 2011, with two earlier observations in 2005 and 2009 (see 
Table [l]). We have analyzed here 20 Swift/XKT observations, 8 RXTE/PCA observations, 
and one CXO/ACIS observation. 



2.1. Swift/XRT data 

Of the 20 Swift/XRT observations listed in Table [TJ four were carried out in the Photon 
Counting (PC) mode and sixteen in the Window Timing (WT) mode which provides much 
better temporal resolution (1.8 ms) at the expense of imaging. We used the HEASOFT0 



1 Version 6.10, 
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analysis tools to reduce and analyze the data. We extracted spectra from the Level 2 event 
data using the standard grade selection of 0-12 and 0-2 for the PC and WT mode data, 
respectively. For the PC mode data, we used an r = 15" circle as the source region and an 
annulus with the same center and inner and outer radii of 30" and 45" as the background 
region. For the WT mode data, we extracted the source spectra using a box centered on the 
CXO location with a length of 30" aligned to the ID image. The background spectra were 
extracted with a similar size box centered far away from the source. We then generated the 
ancillary response files with xrtmkarf for each spectrum and regrouped the source spectra 
with a minimum of 15 counts per bin. The spectral fitting was done in XSPEC 12.6.0. 
Since the source was relatively bright at the onset of the outburst episode, the first XRT 
observation in PC mode (performed during two separate spacecraft orbits) was split into 
two parts to uncover early spectral variations. Three observations in WT mode and three 
observations in PC mode were too short to allow determination of spectral parameters. These 
observations were, therefore, excluded from our spectral analysis. 



2.2. RXTE '/PC A data 

Swift J1834.9— 0846 was observed with RXTE in eight pointings with a total exposure 
time of about 50 ks spanning over 30 days (see Table [l[ ■ The RXTE data were collected with 



the Proportional Counter Array (PC A; Jahoda et al. (1996)) operating with two out of the 



five available proportional counter units in most of the observations. All data were collected 
in the GoodXenon mode, where each photon is time tagged with a minimum time resolution 
of about 1 [is. We used the PCA data primarily for timing analysis as it is not an imaging 
instrument, and the source intensity is relatively dim compared to the bright background 
X-ray emission (e.g., diffuse Galactic ridge emission and bright point sources in the 1° field 
of view of RXTE). However, we extracted the pulse peak spectrum using the longest RXTE 
pointing to investigate the source spectral behavior in a joint PCA and CXO analysis (see 
Section 6.3). 



2.3. CXO data 



We observed Swift J1834.9-0846 on 2011 August 22 with the CXO Advanced CCD 
Imaging Spectrometer (ACIS) operated in the Timed Exposure mode. The target was imaged 
near the aim point on the S3 chip using 1/8 subarray (8' x 1' field of view). The data of an 
archival CXO observation (see Misanovic et al. 2011, for a description) were also analyzed, 



taking into account the different angular resolution and sensitivity. In our analysis we worked 
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with the pipeline-produced Level 2 event files (with standard filtering applied) and utilized 
CIAO 4.3 with CALDB 4.4.5. The spectral fitting was done in XSPEC 12.6.0. 



3. Source Location and Optical Counterpart Search 

We used the wavdetect CIAO tool to determine the point sources in our CXO observa- 
tion. In the vicinity of the Swift/XRT location we find a point source, which we designate 
CXOU J183452. 1-084556, centered at R.A. = 18 h 34 m 52 s . 118, decl. = -08° 45' 56702. 
We also notice the presence of extended emission, up to ~ 15" from the point source, with 
isotropic surface brightness distribution (see Section 6.2.2). The uncertainty of this position 
is dominated by the CXO absolute position uncertainty of 076 (at 90% confidence level) |^] 
The CXO image of the vicinity of Swift J1834.9— 0846 is shown in Figure [TJ 

We compared the CXO image to the archival 2MASS images of the same region of sky. 
We do not detect any near-infrared (NIR) sources within 2" distance from the position of 
CXOU J183452.1 - 084556. We also observed the field of Swift J1834.9-0846 with the Wide 



Field Infrared Camera (WIRC; Wilson et al. 2003) on the 5-m Palomar Hale telescope on 



2011 August 23. WIRC has a field of view of 8'.7 x 8'.7 and a pixel scale of 0.2487 arcseconds 
per pixel. We obtained seven dithered K s band images, consisting of four co-added 30-second 
exposures taken at each dither position. The atmospheric conditions were very good, with 
seeing < 1" and clear skies. The individual frames were reduced in the standard manner 
using IRAF, calibrated, and mosaiced together. The resulting image was astrometrically 
calibrated using 2MASS. The astrometric solution carries a formal la error of 071 for the 
transfer of the 2MASS reference frame to the WIRC image shown in Figure [2j No sources 
are detected within the CXO error circle down to a limiting magnitude of K s ~ 19.5 (at the 
5<T level). The sources designated as SI and S2 on the figure are the ones reported earlier 



by Levan & Tanvir (2011). 



4. Persistent X-ray lightcurve of Swift J1834.9-0846 



Swift J1834.9— 0846 was observed on 20 occasions with Swift after the outburst onset 
(see Table 1). This coverage allows us to construct a lightcurve of the source, which spans 
48 days. In Figure [3j we present the persistent X-ray flux history in the 2— lOkeV range 



as calculated using the power-law (PL) spectral model described in Section 6.1 The X-ray 
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lightcurve of the source indicates a rapid decay in the very early episode (< 1 day), and it is 
consistent with a steady flux decay over the longer term. A PL fit to the temporal decay trend 
(i.e., F oc t~ a ) yields a good fit with a = 0.53 ± 0.03 and a = 0.53 ± 0.07 for the observed 
and unabsorbed fluxes, respectively. Notice that because of the limited spatial resolution, 
the XRT data include both the point source and the surrounding extended emission. As a 
consequence, the decay trend of the point source cannot be unambiguously determined from 
these data. 



5. Timing Analysis 
5.1. RXTE 

Swift J1834.9— 0846 was observed by RXTE on eight occasions with a total exposure 
time of ~ 50 ks, spanning a time baseline of over 30 days (see Table[T]). For our timing analysis 
we used data collected in the 2— lOkeV range. For each observation, we first inspected 
the lightcurve with 0.03125 s time resolution and filtered out the times of short spikes and 
instrumental artifacts. We then converted the event arrival times to that of the Solar System 
Barycenter in Barycentric Dynamical Time (TDB) using the JPL DE200 ephemeris and the 
SW/fc-derived coordinates of the source. 

Next, we employed a Fourier based pulse profile folding technique to determine the spin 
ephemeris of Swift J1834. 9— 0846. We first generated a template pulse profile by folding the 
longest PCA observation (Observation ID: 96434-01-02-00) at the pulse frequency determined 



with a Z\ search (Buccheri et al. 1983). Then, we generated pulse profiles for all PCA 
observations as well as for the CXO pointing, and cross-correlated them with the template 
profile to determine the phase shifts with respect to the template. We obtain the spin 
ephemeris of the source by fitting the phase shifts with a first or higher order polynomial. 
We find that the phase drifts of Swift J1834.9— 0846 are best described with a second order 
polynomial (x 2 =7.3 for 7 degrees of freedom [dof]) that yields a spin period P = 2.4823018(1) 
s and a period derivative P = 7.96(12) x 10~ 12 s s^ 1 (epoch: 55783 MJD). In Figure^ we 
present the drift of the pulse phase with respect to the template and the quadratic trend 
curve (upper panel), and the fit residuals in cycles (lower panel). The measured values of P 
and P correspond to the following spin-down parameters: age r = P/2P = 4.9 kyr, power 
E = 4n 2 IPP- 3 = 2.1 x 10 34 erg s _1 , and magnetic field B = 3.2 x 10 19 (PP) 1 / 2 = 1.4 x 10 14 
G. 



Finally, Figure [5] shows the pulse profiles obtained from several RXTE observations 
folded together using the derived ephemeris. We note the appearance of additional harmonics 
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in the low energy pulse profile of the source (2-5 keV) in RXTE data taken at later times. 



5.2. CXO 

We searched for pulsations in the CXO/ACIS data obtained in the 2011 August 22 obser- 
vation. We used the 733 counts extracted from the r = 1" circle around the CXOU J183452. 1—084556 
position, in the 2-10 keV band (there are only 4 counts below 2keV, likely from the back- 
ground). The time resolution in this observation was 0.44104 s (0.4 s frame time plus 0.04104 s 
charge transfer time). The photon arrival times were transformed to the Solar System 
Barycenter using the CIAO axbary tool. The ACIS observation started at epoch 55795.6489 
MJD and continued for T span = 13.02 ks. 

We calculated the Z\ statistic as a function of trial frequency with a step of 0.35 /iHz 
(which is about 0.057^) and found the maximum Z\ = 467 at v = 0.4028512 Hz±2.0 //HzR 
implying a very high significance of the pulsed signal. We also calculated Z\ for n > 1 but 
did not find a strong contribution of higher harmonics. 



Figure [6] (upper panel) shows the pulse profiles with 5 and 10 phase bins. We used these 
profiles to measure the pulsed fraction)^} p = 85% ± 10%. We estimated the uncertainty of 
the pulsed fraction using Monte Carlo simulations and bootstrapping, also accounting for the 
time resolution and dead time in the 1/8 subarray mode. We also performed randomization 
of the arrival times within the 0.4 s frame time and re-calculated the pulsed fraction, which 
remained within the uncertainty range estimated above. 

The pulsed fraction can also be defined as p = [2(Z% — 2n) /iV] 1 / 2 , where n is the number 
of harmonics that give a significant contribution, and N is the number of count^j In our 
case, p — 1.13 exceeds 100%, which might be due to dead-time effects and the relatively 
large (~ 0.18) ratio of the time resolution to the period. To measure the pulsed fraction 
more accurately, the target should be observed with a better time resolution. 

Figure [6] (lower panel) shows a 20-bin pulse profile averaged over the reference phas^J 



Chang et al. 



2011) 



3 The lcr uncertainty is calculated as 8v — "i 1 / 2 !! x T sp \ n (Z 2 max ) 1 / 2 (see 

4 The pulsed fraction p is defined as the ratio of the number of counts above the minimum level to the 
total number of counts. 

5 The advantage of this definition is the independence of p of phase binning. For the case of purely 
sinusoidal pulsations, p coincides with p (assuming a very large number of bins and low noise), while it is a 
factor of \/2 larger than the RMS measure of variability. 

6 This pulse profile was obtained by averaging 100 pulse profiles (20 bins each) constructed by assigning 
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We also produced a pulse profile for the surrounding extended emission but did not find a 
statistically significant pulsed signal. 



6. Spectral Analysis 
6.1. Swift 

We have fitted all XRT spectra (2 — 10 keV) jointly with two continuum models: a power 
law (PL) and a single blackbody (BB), both with interstellar absorption. In the first case, we 
find that the photon index remains the same within the uncertainties; therefore, we forced all 
observations in our joined fit to have the same varying photon index, while the normalizations 
were allowed to vary individually. We obtained a good fit (xt = 1-01 f° r 62 dof) with the best 
model parameters iVn = 10.5^'g x 10 22 cm -2 , and photon index T = 3.2 ±0.4. The absorbed 
BB model resulted in temperatures that also remained consistent within their uncertainties; 
we then linked the temperatures and allowed the normalizations to vary. We again obtained 
a good fit (xl = 1-04 for 62 dof) with N H = 4.4l{;^ x 10 22 cm" 2 and kT = 1.1 ± 0.1 keV. 
The temperature is higher than those measured in most other magnetars (typically around 



0.5 keV; Woods & Thompson 2006) 



6.2. CXO 

6.2.1. CXOUJ183452. 1-084556 

We collected a total of 733 counts (2-10 keV) from a circular region of r = 1" centered 
at CXOU J183452. 1-084556; the background contribution is expected to be only 0.25 counts 
(background was measured an 20" < r < 33" annulus). We then grouped the source spectrum 
requiring a minimum of 15 counts per spectral bin. The resulting spectrum is shown in 
Figure [7] (black error bars). The source pileup is negligibly small (< 1%), as the total 
source count rate of 0.057 counts s -1 corresponds to 0.025 counts per frame. The point 
source spectrum can be fitted equally well with both the absorbed PL and BB models (see 
Table ||. The observed (absorbed) flux (2 - lOkeV) is F point = (3.0 ± 0.1) x 10~ 12 erg s" 1 
cm -2 (corrected for the finite extraction aperture and 9% deadtime). Table [2] contains the 
values of the iVn and photon index for the best-fit absorbed PL model, and the Nr and 
temperature (kT) of the BB model. From the BB fit we estimate the emitting area radius 



different phases to the first count and folding with the SGR period. 
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to be 0.26 km (assuming that the source is at the same distance of 4kpc as the SNR W41; 



Tian et al. 2007). The corresponding unabsorbed PL and BB fluxes (2 — lOkeV) are 1.6 



+0.6 
-0.4 



and (5.8 ± 0.6) x 10 erg s , respectively (see also Table |2J). 

Guided by the pulse profiles shown in Figure |6j we have extracted spectra from two 
different phase intervals: 0.15-0.35 (peak) (indicated with the shaded region in Figure [6 
(upper panel)) and the rest (off-peak). These spectra are shown in Figure |8j we again used 
the PL and BB models in each case, fixing the Nr at the best-fit value of the phase-integrated 
spectrum. The peak and off-peak spectral parameters (r or kT; see Table [3]) are consistent 
within their uncertainties. 



6.2.2. Halo 

We collected 314 counts (2-10 keV) from the 2" < r < 10" annulus (hereafter "halo"), 
centered at CXOU J183452. 1—084556, where the contribution of the point source is expected 
to be small (< 10%). We subtracted the background (estimated from a much larger region 
away from the source) and obtained a net total of ~ 300 counts. To separate the halo from the 
point source, we simulated a point spread function (PSF) using MAPuXQ The comparison 
of the data with the PSF simulation (Figure [9]) shows a good agreement within a small 
aperture (approximately up to 1" radius), while the extended emission dominates at larger 
radii. Based on our simulation, we estimate that ~ 33 photons come from the point source, 
after taking into account the extended PSF wings. The final halo spectrum was also binned 
requiring a minimum of 15 counts per bin. The best-fit PL slope is approximately the same 
as that of the point source spectrum, while the best-fit Nn is a factor of two lower (see Table 
2]). The A^H-r confidence contours for the halo spectrum, together with those for the point 
source spectrum, are shown in Figure [lOj The absorbed and unabsorbed fluxes (2-10 keV) 
of the halo emission are F halo = (4.7 ± 0.2) x 10~ 13 and (1.6 ± 0.4) x 10~ 12 erg s _1 cm -2 , 
respectively 

The extended emission is well described by the same dust halo model as the one used 



by Misanovic et al. (2011) for another nearby source, CXOU J183434. 9— 084443, according 
to which most of the dust must be located relatively close to the source (within 1/4 of 
the distance). At least part of this dust could be associated with the molecular cloud that 



appears to be interacting with W41 (Leahy & Tian 2008), in agreement with the very large 



absorption column that we find in our spectral analysis (see Table 2). 



http : / / space . mit . edu/ cxc/marx/ 
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6.2.3. Pre- outburst CXO/ACIS data taken on June 2009 
We analyzed the 2009 CXO observation covering the Swift J1834.9— 0846 field and found 



zero photons within the error circle (r = ff'.Q) of CXOU 183452.1—084556 (see Figure 11). 
The off-axis angle of ~ 4'.6 during that observation is, however, large enough for the angular 
resolution to be substantially degraded compared to on-axis. Hence, to estimate the 2009 
upper limit on the source flux, we used a larger radius, r = 2", which would contain about 
50% of the flux of a point source at this off-axis angle. We found 5 and 4 photons in the 
0.5-10 and 2-10 keV bands, respectively. The mean local background surface brightness is 
0.24 ± 0.02 and 0.18 ± 0.01 counts arcsec" 2 in 0.5-10 and 2-10 keV, respectively. Thus, 
within the r = 2" extraction aperture we would expect to detect about two counts from the 
background. This translates into an upper limit of 0.15 counts ks _1 in 2 — lOkeV, which 
corresponds to an absorbed flux limit of Fp 0int < (2 — 4) x 10~ 15 erg s _1 cm -2 . 

Although we do not detect a point source, we notice extended emission on larger scales 
around the position of Swift J1834.9— 0846. We demonstrate this by plotting the radial profile 



of the surface brightness (see Figure 11 , top panel). The ten annuli used to extract the radial 
profile are centered at the position of Swift J1834.9— 0846 , while the background is measured 
from ten circular (r = 20") regions surrounding the source. One can see from the figure inset 
that most of the excess over the background is within r < 12" and it corresponds to a de- 
tection significance of ~ 5.1cr. However, there is also marginal (~ 3a) evidence for extended 



emission at larger scales (between r = 12" and 30", see Figure 11, top). There might be 
an even more extended (primarily toward southwest from the Swift J1834. 9— 0846 ), fainter 
asymmetric emission, but its significance can only be established with deeper CXO/ACIS 
observations. 

We also found evidence for extended emission in the 2005 XMM-Newton data (see also 



Mukherjee et al. 2009). The extent, location and significance of the extended emission in 
the EPIC/MOS images (which are not affected by the chip gaps and have a low enough 
background) are similar to those measured from the 2009 CXO/ACIS images (see above). 
The previously reported large-scale extended emission west-northwest of the SGR (i.e., in 
the direction toward CXOU183434. 9-084443) could be mainly due to the point sources that 
are smeared out in the EPIC images because of the coarse angular resolution of XMM- 
Newton. The two brightest point sources are clearly resolved in the sharper CXO/ACIS 



images (cf. corresponding panels in Fig. 11). We did not attempt to extract the spectra from 
the 2005 XMM-Xewton/EPIC data because the background is much higher and the angular 
resolution is worse than the one of the CXO observation. No point source is detected in the 
XMM-Newton images at the position of CXOU J183452.1-084556. We have not estimated 
an XMM-Newton upper limit on the point source flux, as it would be less restrictive than 



-li- 



the one derived using the 2009 CXO/ACIS data. 



6.3. Joint fits of CXO and RXTE data 

Since the RXTE/PCA is not an imaging instrument, we could not spatially separate 
the halo and the point source or even subtract a background measured independently from 
an offset region. However, since the instrument has a broader spectral range than the 
Swift/XKT and CXO/ACIS, potentially providing valuable source information above 10 keV, 
we used the latter data to calibrate our PCA spectrum of the longest (9.7 ks) pointing of 
2011 August 9 (see Table [I]). For the RXTE data, we accumulated the spectrum at the pulse 
minimum (which contains background, halo, and any unpulsed point source contributions) 
as the background and subtracted it from the source spectrum integrated over the remaining 
phases (see shaded regions in Figure [5J. The resulting pulsed emission spectrum was then 
rebinned to have at least 50 counts per spectral bin after the background subtraction. 

We performed a joined spectral fit of the RXTE (2-50 keV) and the CXO/ACIS (2- 
lOkeV) data (see Figure [7]). We found that the best fits are obtained when the RXTE flux 
is scaled down by a factor rj = 0.6. The resulting PL best-fit parameters are very close to 
those of the CXO/ACIS fits but somewhat better constrained (see Table [2]). A single BB 
fit is disfavored by systematic residuals at energies > 8keV (see Figure [7j bottom panel). 
The introduced scaling of the RXTE flux can be interpreted as due to two reasons: (1) the 
source was brighter at the time of the RXTE observation, and (2) the true background is 
lower than that estimated from the pulse minimum (see above). 



7. Discussion 



7.1. Swift J1834.9-0846 



Swift J1834.9— 0846 has one of the shortest periods among magnetar^] and one of the 
highest pulsed fractions of the persistent X-ray emission, similar to that of IE 1048.1—5937 on 



2000 December 28 (Tiengo et al. 2002 ). In its timing and spectral properties, Swift J1834.9— 0846 



strongly resembles the recently discovered SGR J1833— 0832, which has a period of 7.6 s and 
a magnetic field of 1.8 x 10 14 G ( Gogii§ et al.|2010 Esposito et al.|20il~ ). In particular, similar 



*See the McGill AXP/SGR catalog: 



http : //www . physics . mcgill . ca/$\sim$pulsar/magnetar/main . 



html 
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to SGR J1833— 0832 (and unlike most other SGRs with good quality spectra), the 0.5—10 keV 
spectrum of Swift J1834.9— 0846 can be fitted with a single BB model, whose temperature, 
kT ~ 1.1 keV, is the same as that of SGR J1833-0832. The Swift J1834.9-0846 BB radius, 
R = 0.26 km, is a factor of 3 smaller than in SGR J1833— 0832, which, however, may not 
be a significant difference given the poorly known distances. Another similarity between 
Swift J1834. 9— 0846 and SGR J1833— 0832 is the lack of obvious spectral shape evolution 
with rotational phase. The phase-resolved spectra (see Figure |sj) differ only in normaliza- 
tion, and the differences in other model parameters are not statistically significant. Despite 
these similarities, the post-burst flux decay trend is markedly different for the two SGRs. 
The unabsorbed flux of Swift J1834. 9— 0846 decreased as cx £- a53±ao7 from day 2 after the 
burst (Figure [3]), while the flux of SGR J1833— 0832 remained constant for nearly 20 days 
before the onset of decline. We note, however, that this early constancy of the flux in 
SGR J1833— 0832 is unusual; the enhanced persistent X-ray flux of magnetars following an 
outburst usually declines as a power law with an index similar to that of Swift J1834.9— 0846. 

Similar BB temperatures and radii were also found for SGR 0418+5729 (P = 9.1 s, 
B < 7.5 x 10 12 G; |Esposito et~aT1|2010[ ) from the Swift XRT data taken within ~ 10 days 
after the outburst. Also, SGR 0418+5729 exhibited a cx t~ 0,3 decay during the first 19 days 
and a much steeper, oc t~ 1 ' 2 , decay thereafter. 

It is tempting to interpret the small emitting area of Swift J1834.9— 0846 (similar to 
those of SGR J1833— 0832 and SGR 0418+5729) as a hot spot on the neutron star surface. 
We should note that it would be very difficult to obtain such a high pulsed fraction even for 
a very small hot spot emitting (nearly isotropic) BB radiation because the pulsations would 



be washed out by the light bending in the neutron star gravitational field (see, e.g., Zavlin et 



al. 1995). If, however, we take into account that the angular distribution of radiation from 



a neutron star atmosphere has a narrow peak along the magnetic field direction (Pavlov et 



al. 1994), such a high pulsed fraction can indeed be explained assuming that the observed 
radiation emerges from a small hot spot near the magnetic pole of the neutron star. The fact 
that the light curve with such a high pulsed fraction shows only one peak per period suggests 
that the magnetic field configuration is substantially different from a centered dipole (e.g., 
it could be a strongly decentered dipole, in which case the magnetic fields and the temper- 
atures are substantially different at the two poles). We caution here that the BB model 
provides only an empirical description of the spectral shape. It can be used for comparison 
of different sources, but it may be significantly different from the actual spectrum emitted 
from a neutron star atmosphere ( Pavlov et al.||199~5 ) and possibly modified by the resonance 



Compton scattering in the neutron star magnetosphere (Nobili et al. 2008). Comparing the 
BB fit parameters of the three SGRs (J1834.9-0846, J1833-0832, and J0418+5729), we can 
conclude that they do not depend on the SGR period (in the range of 2-10 s), nor on the 
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strength of the spin-down magnetic field (in the range of (0.1-2) x 10 14 G). 

The nondetection of the Swift J1834.9— 0846 in the pre-outburst CXO data shows that 
the SGR flux can vary by at least a factor of ~ 10 3 between the presumably truly quiescent 
level in the low state and an elevated level that has persisted, with a slow decay, for at least 
6 weeks after the outburst. This suggests that there is a large number of SGRs in a quiescent 
state undetectable at the current level of sensitivity of X-ray observatories. 



7.2. Extended emission 

At first glance, the extended emission around Swift J1834.9— 0846, detected by CXO in 
2011, looks rather unusual. Its radial distribution is consistent with that of a dust scattering 
halo (see Figure [9]), but its spectrum shows some peculiarities. In particular, the best-fit 
hydrogen column density of the halo is a factor of 2 lower than that of the central source, while 
their spectral slopes are similar within statistics, instead of being steeper by Ar = 1-2, as 
expected for the model halo spectrum (see, e.g., Misanovic et al.|[2011 ). A likely explanation 



can be derived from the breadth of the Nu-T confidence contours and the strong correlation 
of these parameters. Indeed, Figure 10 shows that an intermediate Nn ~ 1.5 x 10 23 cm -2 



corresponds to the 90% point source and halo confidence contours, and the best-fit photon 
indices at such TVh are r ~ 3 and ~ 5 for the point source and halo, respectively. Thus, 
we believe that the dust scattering halo is the most plausible interpretation of the extended 
emission around Swift J1834.9— 0846. 

We should also note that a fainter extended emission was seen around the magnetar 
positition in the archival CXO data from 2009, in which no point source was detected. 
Although the best-fit and T for this pre-outburst emission are substantially smaller than 
those in the post-outburst data (see Table [2]), the large uncertainties of these parameters 
make them consistent with the corresponding parameters measured for the halo in 2011. The 
existence of a halo in the archival data may indicate that Swift J1834. 9— 0846 experienced 
an outburst not long before the 2009 June 7 observation. 

Although it seems certain that most of the extended emission is the dust-scattered 
emission from the magnetar, we cannot exclude the possibility that it may contain some 
kind of a pulsar wind nebula (PWN), due to synchrotron radiation from relativistic elec- 
trons/positrons accelerated in the neutron star magnetosphere and shocked in the ambient 
medium. We know that in the case of rotation-powered pulsars, a typical X-ray PWN lumi- 



nosity is about 10 E 1 albeit with a large scatter (Kargaltsev & Pavlov 2008). If the same 



relationship is valid for magnetars, we would expect L pwn ~ 10 30 erg s 1 , which would be 
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undectable at the presumed distance of 4kpc. It might happen, however, that a "magnetar 
wind nebula" is more efficient than one created by a rotation-powered pulsar, in which case 
we would expect a detectable contribution. To separate it from the dust scattering halo, one 
should analyze several data sets obtained at different times after the outburst. We expect 
that the halo component flux would be changing in proportion to the point source flux (with 
a time lag), while the PWN component would remain constant. 



7.3. Relation to SNRW41 and HESS J1834-087 



The distance to Swift J1834.9— 0846 still remains an open issue. As this source is lo- 
cated within W41, association with this SNR is certainly plausible (other SGRs were found 
near SNR centers, |Woods fc Thompson 2006), but it has not been firmly proven. Similar 
extreme absorption (JVh = 3 x 10 23 cm -2 ; Misanovic et al. 2011) has been measured for the 
neighboring CXOU J183434. 9— 084443, indicating that such an absorption is not a unique 
feature of Swift J1834.9— 0846 and hence not intrinsic to it. However, the distance to (and 
the origin of) CXOU J183434. 9— 084443 are also unknown. It could be a pulsar associated 



with W41 or a background AGN located much farther. Leahy & Tian (2008) presented evi- 



dence for molecular clouds near W41, which are likely interacting with the SNR. The large 
absorbing column could be attributed to those clouds. At this point we can only conclude 
with certainty that Swift J1834.9— 0846 is at the distance of 4kpc or farther (the line of sight 
in that direction intersects several spiral arms). To better constrain the distance, the method 



of Durant & van Kerkwijk (2006) could be used; however, it requires grating observations, 



which are only feasible when the source is in the bright state. We note, however, that the 



conclusion by Durant & van Kerkwijk (2006), that all AXPs have more or less standard lu- 



minosity of 1.3 x 10 erg s , cannot hold for quiescent SGRs because otherwise they would 
have been easily seen even at the most extreme distance of 20 kpc. 

The field surrounding Swift J1834. 9— 0846 is rich with high-energy sources (see Figure 
13). This magnetar is located at the heart of SNRW41 and nearly at the center of the 
extended TeV source HESS J1834— 087, which itself is confined to the SNR interior (see 



Figure 13 ). In addition, there is a somewhat offset Fermi source, 2FGL J1834.3— 0848, located 



nearby (see Figure 13 ). Since the extent of HESS J1834— 087 is significantly smaller than that 
of the SNR, the TeV emission cannot be coming from the SNR shell as it does in some other 



cases (Bochow et al. 2011). The only other plausible explanation is that the TeV emission 



is powered by relativistic electrons injected by the compact object formed after the SNR 
explosion. There are currently several candidates for such an object. Firstly, a few SGRs 



are known to be associated with shell-type SNRs (Hurley 2000), and the central location 



- 15 - 



of Swift J1834.9— 0846 certainly supports such a hypothesis. On the other hand, there is no 
firm evidence so far that SGRs can produce copious amounts of relativistic particles similar 
to young rotation-powered pulsars. While there is a convincing evidence that pulsars can 
power relic PWNe emitting TeV 7-rays, such evidence is currently lacking for magnetars. 

Among other sources possibly related to HESS J1834-087 and W41 are CXOU J183434.9- 



084443 (a PWN candidate discussed in detail by |Misanovic et al.||20lT| ) and the 2XMM 
J183417.2— 084901, which is located right at the center of the Fermi error circle (see Fig- 
12, top panels). Further longer observations of this region are required to understand 



ure 



the connection between the sources observed in different energy domains. 
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Table I. X-ray observations of Swift J1834.9-0846. 



Date 


ObsID 


Observatory/Detector (Mode) 


Exposure, ks 


Time resolution, s 


2005 Sept 18 


0302560301 


XMM-Ncwton EPIC 


18.6 


0.072 


2009 Jun 7 


10126 


CXO ACIS-S 


46.5 


3.2 


2011 Aug 7 


00458907000 


Swift/XRT (PC) 


1.54 


2.5 


2011 Aug 7 a 


00458907001 


Swift/XRT (WT) 


0.096 


1.8 x 10~ 3 


2011 Aug 8 a 


00458907002 


Swift/XRT (WT) 


0.129 


1.8 x 10~ 3 


2011 Aug 8 


00458907003 


Swift/XRT (WT) 


1.65 


1.8 x 10~ 3 


2011 Aug 8 


00458907004 


Swift/XRT (WT) 


0.958 


1.8 x 10~ 3 


2011 Aug 9 


00458907006 


Swift/XRT (WT) 


2.67 


1.8 x 10~ 3 


2011 Aug 9 


96434-01-01-00 


RXTE PCA 


3.40 


9 x 10~ 7 


2011 Aug 9 


96434-01-02-00 


RXTE PCA 


9.66 


9 x 10~ 7 


2011 Aug 12 


r\r\ a tr onriTnnv 

U0458907007 


bwtft/XRL (Wl) 


5.67 


1.8 X 1U 


2011 Aug 14 


00458907008 


Swift/XRT (WT) 


5.39 


1.8 X 10~ 3 


2011 Aug 14 


96434-01-03-00 


RXTE PCA 


6.78 


9 x 10~ 7 


2011 Aug 18 


96434-01-03-01 


RXTE PCA 


6.75 


9 x 10~ 7 


2011 Aug 18 


00458907009 


Swift/XRT (WT) 


5.73 


1.8 x 10~ 3 


2011 Aug 21 


00458907010 


Swift/XRT (WT) 


2.49 


1.8 x 10~ 3 


2011 Aug 22 


14329 


CXO ACIS-S 


13.0 


0.44104 


2011 Aug 24 


96434-01-04-00 


RXTE PCA 


6.60 


9 x 10- 7 


2011 Aug 24 a 


00458907011 


Swift/XRT (WT) 


0.94 


1.8 x 10~ 3 


2011 Aug 27 


00458907012 


Swift/XRT (WT) 


1.95 


1.8 x 10" 3 


2011 Aug 29 


96434-01-05-00 


RXTE PCA 


6.05 


9 x 10~ 7 


2011 Aug 30 


00458907013 


Swift/XRT (WT) 


2.16 


1.8 x 10~ 3 


2011 Sep 2 


96434-01-06-00 


RXTE PCA 


5.12 


9 x 10" 7 


2011 Sep 2 a 


00458907014 


Swift/XRT (PC) 


2.06 


2.5 


2011 Sep 5 a 


00458907015 


Swift/XRT (PC) 


1.72 


2.5 


2011 Sep 8 


96434-01-06-01 


RXTE PCA 


5.52 


9 x 10- 7 


2011 Sep 10 a 


00458907016 


Sun/S/XRT (PC) 


2.01 


2.5 


2011 Sep 15 


00032097001 


Swi/i/XRT (WT) 


9.09 


1.8 x 10~ 3 


2011 Sep 18 


00032097002 


Swi/tyXRT (WT) 


10.45 


1.8 x 10~ 3 


2011 Sep 21 


00032097003 


Swift/XRT (WT) 


7.44 


1.8 x 10~ 3 


2011 Sep 24 


00032097004 


Sm/VXRT (WT) 


8.10 


1.8 x 10~ 3 



Note. — Log of all observations used in our analysis. 
a Excluded from the spectral analysis. 
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Table 2: Spectral fit parameters of Swift J1834.9— 0846 (point source) and its extended emis- 



sion (halo). 

Model ' iy H ,22 ~ Af a orR h T c or kT c x^/d.o.f. d L x c 





CXO, PL (point) 


20.3l|£ 


4.81H 


q pr+0.5 


0.72/31 


3.1 




CXO, BB (point) 


12.01};* 


0.2618$ 


1.1 ±0.1 


0.75/31 


0.33 




CXO, PL (halo) 


10.01^ 




3.7 ±0.6 


0.86/17 


0.31 




CXO, PL (point) 


15 (fixed) 


0.8318;" 


2.7 ±0.2 


0.86/32 


1.7 




CXO, PL (halo) 


15 (fixed) 


4.8±?:| 


5.0 ±0.3 


1.05/18 


0.67 


cxo, 


PL (pre-outburst, halo) 


4.0l?i 


5.5±|;g x 10~ 4 




0.72/7 


~0.01 


cxo, 


PL (pre-outburst, halo) 


15 (fixed) 


4.7+f 2 x 10~ 2 


3.5 ±0.1 


0.80/8 


-0.03 


cxo 


(point) and RXTE, PL 


20.7ll;g 


4.81^ 


q fi + 0.4 


0.77/41 


2.8 


cxo 


(point) and RXTE, BB 


lo.olrl 


0.181^° 


i i pr+0.09 
i - io -0.08 


0.88/41 


0.18 



Note. — The uncertainties are given at 68% confidence level for a single interesting parameter. 
a Spectral flux in units of 10 -2 photons cm~ 2 s _1 keV^ 1 at 1 keV. 
b BB radius, in units of km 2 . 
c Photon index or BB temperature in keV. 
c Rcduccd x 2 an d the number of degrees of freedom.. 

e Unabsorbed PL luminosity in the 2—10 keV band or bolometric BB luminosity (irR 2 aT 4 ), in units of 10" 
ergs s . 



Table 3: Power-law fits to the phase resolved ACIS spectra of Swift J183 4.9— 0846. 



Phases 


^H,22 


A/" a 


r b 


tf/d.o.f. c 


peak 


20.3 


1.5 ±0.2 


3.8 ±0.3 


0.94/40 


off-peak 


20.3 


0.30 ±0.04 


3.5 ±0.2 


0.94/40 



Note. — Nh was held fixed during the fit. The uncertainties are given at 68% confidence level for a single 
interesting parameter. 
"Spectral flux in units of 1CT 2 photons cm~ 2 s _1 keV -1 at 1 keV. 
b Photon index. 

c Reduced x 2 and the number of degrees of freedom. 
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Fig. 1. — Image of CXOU 183452.1-084556 and surrounding emission (0.7—10 keV) obtained 
with ACIS-S3 on 2011 August 22. The radii of the inner and outer circles are 2" and 12", 
respectively. 
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Fig. 2. — Palomar/WIRC i^ s -band image showing the r = 0.6" CXO error circle for 
CXOU 183452.1-084556. The sources designated as SI and S2 are the ones reported by 
Levan & Tanvir (2011). 
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Fig. 3.— Persistent X-ray lightcurve (2 - lOkeV) of Swift J1834.9-0846/CXOU J183452.1- 
084556 obtained from 48 days monitoring of the source with Swift /XRT. The dashed line 
shows the best-fit power-law temporal decay model (oc t _a53 ). 
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Fig. 4. — Top panel: Plot of phase shifts for each RXTE observation of Swift J1834.9— 0846. 
The solid line is a quadratic trend that fits the time evolution of the phase shifts. Bottom 
panel: Residuals of the fit. 



Aug 9 - Aug 14 



Aug 18 - Sep 8 




Fig. 5. — Dependence of the RXTE pulse profiles on time and energy. The shaded area in 
the top left panel corresponds to the phase interval used for spectral analysis (see Section 
6i3b. 
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Fig. 6. — Top panel: CXO pulse profiles (2 — lOkeV) with 5 and 10 phase bins. The shaded 
regions indicate the peak interval (phases 0.15-0.35) used for phase-resolved spectroscopy. 
Bottom panel: CXO pulse profiles (2 — 10 keV) with 10 phase bins, averaged over the reference 
phase. 
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Fig. 7.- CXO/ACIS (black/bottom) and RXTE/PCA (red/top) spectra of 
CXOU 183452.1-084556 jointly fitted with the PL and BB models (top and bottom panels, 
respectively). 
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Fig. 8.— PL fits of the CXO spectra of CXOU 183452.1-084556 in the pulse maximum 
(black/bottom) and pulse minimum (red/top). 
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Fig. 9. — Radial profile in 2 — 10 keV from the 13 ks CXO/ACIS observation (histogram) 
shown together with the simulated PSF (dashed line with error bars), a dust halo model 
(dashed line) from Misanovic et al. (2011), and a background (horizontal dash-dotted line) 
measured from the current observation. The solid line shows the sum of all three components. 
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Fig. 10. — Confidence contours (68% and 90%) in the Nr-V plane for the PL fit to the 
halo (blue), point source (green), and the pre-outburst diffuse emission (black) spectra. The 
contours are obtained with the PL normalization fitted at each point of the grid. The best-fit 
parameter values are shown by crosses. 



Fig. 11 — Archival 2-8 keV ACIS-S3 images of the CXOU J183452. 1-084556 field obtained 
on 2009 June 7. Upper panel: Binned (pixel size C/98) and smoothed (with a 2'.'9 Gaussian 
kernel) image of the field. The annuli, centered at the SGR position, are used to extract the 
radial profile of the surface brightness distribution (shown in the inset). The background 
(solid horizontal line in the inset) and its uncertainty (dotted horizontal lines) are measured 
from the ten r = 20" circular regions. Lower panel: Zoomed-in image of the field around 
the position of CXOU J183452. 1-084556 at the native ACIS-S3 binning with no smoothing 
applied. The circles with radii 0'.'6 and 2" were used to estimate the detection upper limit 
(see Section 6.2.3). 
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Fig. 12. — Mosaic of images (2-10 keV) of the central region of the W41 SNR, spanning 
six years. The left, middle, and right columns correspond to the data obtained with XMM- 
Newton/EPlC (2005 Sept 18), CXO/ ACIS-S (2009 June 7), and Swift /XRT (2011 August 
7), respectively. The lower panels are zoomed in the Swift J1834. 9— 0846 position. 



Fig. 13. — Multiwavelength emission from the W41 region. The left panel shows the CXO 
ACIS image (0.3-8 keV; red) and the VLA 20 cm image (blue) from the MAGPIS database 
(http://third.ucllnl.org/gps/). In the right panel the 20 cm VLA image is shown in blue, 
the Spitzer IRAC 8 jum image in green, and the Spitzer MIPS 24 fim image in red. The 
r = 5'.4 circle shows the extent of HESS J1834-087 ( Aharonian et al.]|2005 ); the ellipse shows 
the position of 2FGL J1834.3-0848 at the 95% confidence level (lAbdo et al.||2011b. 



